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STATISTICAL DISTRIBUTION OF THE BREAKING STRENGTH OE A 


BUNDLE OF CLASSICAL FIBERS 

by Burt M. Rosenbaum 
Lewis Research Center 


SUMMARY 

It is shown that, for ideal bundles having a given number of fibers that 
will break under a load per fiber Zp and a given number of fibers that will 
not break under a load Zg (>^p), the average number of unbroken fibers in the 
bundle at loads between Zp and Zg depends linearly on the probability that 
an individual fiber picked at random will support a load Z. Bundles are 
grouped according to the number of fibers in the bundle that can support that 
load per fiber at which the average total load supported by the bundle is a 
maximum. An approximation to the maximum of the average load supported by the 
bundle is arrived at for each group, and this maximum averaged over the groups 
constitutes a lower bound for the average breaking strength of the bundle. The 
results agree with those of Pierce and Daniels. 


INTRODUCTION 

A classical fiber is one that will support a load less than its breaking 
strength indefinitely without stretching or breaking but will break immediately 
under any load equal to or in excess of its breaking strength. The fundamental 
papers dealing with the breaking strengths of an ideal bundle (no friction or 
twist) of equal-length classical fibers are those by Pierce (ref. l) and Daniels 
(ref. 2). Pierce was the first to obtain an approximation to the average break- 
ing strength of a bundle of N such fibers. Daniels rigorously showed by 
sophisticated mathematical methods that the expression of Pierce is the correct 
asymptotic limit of the breaking strength as N gets large. He also found the 
asymptotic limit of the variance of the breaking strength for large N. In 
this report, the method of reasoning employed leads to the establishment of the 
aforementioned results from a viewpoint that is based on physical rather than 
mathematical insight. 


SYMBOLS 

B( Z ) probability that randomly selected fiber will break at load less than Z 

b(Z) frequency function for breaking strength of individual fibers randomly 
selected 

E expectation value 



G( Z ) ZR(Z) 


G n derivative of G( Z ) with respect to Z evaluated at Z 

H n (y) Hermite polynomial of y 

k dummy index 

L total load on bundle of fibers 

L( Z ) total load at Z averaged over population of bundles 
Z load per fiber 

Z load per fiber at which average total load L is maximum 

N number of fibers in bundle 

n number of unbroken fibers in bundle 

n number of unbroken fibers in bundle at load Z 

0 order 

P probability 

R(Z) probability that randomly selected fiber will support load up to Z 

S average breaking strength of bundle 

y variable (NR(z) - n) j ^NR( Z)B(Z) 

S£ (NR(Z) - n)/(N - n) 

9 variable of integration 

A£ value of Z that maximizesjaverage load supported by bundles having n 

unbroken fibers at load Z 

l-i^y third moment about origin for variable y 

a standard deviation 

CONCEPT OF PATH 

Let the distribution of breaking loads of a classical fiber be given by 
the frequency function b(Z), where b(Z)dZ is the probability that a randomly 
selected individual fiber will have a breaking strength lying in dZ at Z . 
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The cumulative distribution function 

b(i) = f b(e)de 

J o 

is the probability that a randomly selected fiber will break at some load less 
than Z. Let R(Z)=1-B(Z) denote the probability that a fiber will remain 
unbroken up to load Z . 

When a load L is hung on an ideal bundle of N classical fibers , the 
load distributes itself evenly over the unbroken fibers so that each of the n 
unbroken fibers supports a load Z = L/n. The breaking strength of a given 
fiber of the bundle does not depend on the action of the other fibers of the 
bundle and, hence , the breaking strength of the bundle may be determined from 
the strengths of the individual fibers. In particular, suppose one considers 
the bundle as separated into its independent individual fibers with separate 
loads Z of equal magnitude on each fiber. Let this load Z be allowed to 
increase monotonically from zero to the breaking strength of the strongest 
fiber. Then, at any point in the process where the individual fibers support a 
load Z, the total load L( Z ) supported by the bundle is the product of the 
number of unbroken fibers n( Z ) and the load Z acting on each unbroken fiber. 
Thus, the breaking strength of the bundle is the maximum of the total load L. 
After this maximum L has been reached, further increases in Z cause a reduc- 
tion in the total load L( Z ) = Zn(Z). (if the fibers were not separated in this 
manner, and if breaking strength were determined by increasing L until break- 
age occurred, no decreases in total load L would appear, and the bundle would 
immediately break after the breaking strength of the bundle had been reached. ) 

The situation may be depicted by means of figure 1 (p. 4). In fig- 
ure 1(a) , a typical plot of Z against R( Z) is given. (Note that in the fig- 
ures, R( Z) decreases from unity to zero as one proceeds from left to right along 
the horizontal axis.) It can be seen that as Z increases from zero to high 
values, R( Z) decreases from unity to zero or, if R(z) is to be considered the 
independent variable, Z increases from zero to high values as R(z) decreases 
from unity to zero. The path followed by a bundle of N fibers averaged over 
the population of bundles is shown as the straight line in figure 1(b), where 
the expected number of surviving fibers is given by E(n) = NR(Z). The path 
taken by a typical bundle is shown in figure 1(c), where n for that bundle 
is plotted against R(Z). For large N at least, the actual path followed by 
a particular bundle would be expected to lie fairly close to the average path 
as given by the E(n) line. 


RESULTS OF PIERCE AND DANIELS 

The average load supported by the bundle as a function of Z is 

L(l) = l E(n) = WZR(Z) (l) 

The maximum average load supported by the bundle is 

W = L(Z) = NZR(2) (2) 
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where Z 



satisfies the equation 



(c) Number of unbroken fibers. 

Figure 1. - Plots with R (l) as independent variable. 


dh( Z ) 
dZ 


= 0 


Z=Z 


R(Z) = Zb( Z ) 


( 3 ) 


Equations ( 2 ) and (3) represent the 
approximation to the breaking strength 
of the bundle as given by Pierce, 
Daniels showed that this result is the 
correct asymptotic limit for the ex- 
pectation value of the breaking 
strength for large N. He also found 
that the asymptotic limit for the 
standard deviation of the strength for 
large N is given by 


c = Z 


|/kr(z)b( 


Z) 


(4) 


Daniels makes the pertinent remark that 
equation (4) would follow if one could 
assume that the breaking strength is 
only dependent on the number n^= n(Z) 
of fibers surviving at a load Z where 
n is, of course, distributed according 
to' the simple binomial law. He adds, 
however, that there appears to be no 
"a priori justification” for this 
assumption. 

The average breaking strength S 
of a bundle of N fibers is the maxi- 
mum load that the bundle supports 
averaged over the population of bundles . 
For a particular bundle, this maximum 
load^will not, in general, occur at 
Z = Z, and hence, it will be larger 
than L( Z ) , the value of L ^at Z . On 
the other hand I^ax is L( Z ) averaged 
over the population of bundles, and 
therefore, constitutes a lower 


bound for S. If the maximum load for every_possible path (in the_sense of 
fig. 1(c)) occurred at Z = Z, then Lmax = L(Z) would indeed be S exactly 
It can be said that the approximation to S as given by equation (2) treats 
every path as though its maximum load occurred at Z = ?. 


AVERAGE PATH BETWEEN TWO DESIGNATED POINTS 
It shall first be shown that, for values of Z lying between Zq and 
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the average path taken by bundles which are constrained to go through the 
points R(Zp), np and R(Z 2 ),n 2 (i.e., bundles which have np unbroken 
fibers at a load Zp and n2 unbroken fibers at a load Z 2 ) is merely the 
straight line connecting these two points on a plot of n against R( Z ) . The 
proof is straightforward and goes as follows. 

The probability that a bundle originally consisting of N fibers takes 
a path such that it has np survivors at Zp and n2 survivors at Z 2 
where IT > np > n 2 > 0 and 0 < Zp < Z 2 is 


P(np at Zp, n 2 at Z 2 ) 


ITl 

(N - np) 1 (np - n 2 )Jn 2 i 


[B(Zp)] W - n l[R(Zp) - R(Z 2 )] n l- U 2[R(Z 2 )] n 2 (5) 

The probability that a bundle originally consisting of N fibers takes a path 
that passes through the three points np at Zp, n at Z, and ng at Z 2 , 
where ^>n^>n>n 2 >0 an & 0 £ ^1 £ ^ — Z 2 , 


P( ni at l 1; n at I , n 2 at l 2 ) - ( M . ni ) : (^ -"n) i ' (n - ' n 2 ) !n 2 ! 

[b(1 1 )] N - D 1[e (I i ) - R (! )] n l- n [R(D - R(l 2 )] n - n 2[R» 2 )P (6) 


The conditional probability that a bundle which passes through np at Zp and 
n 2 at Z 2 has n unbroken fibers at Z is given by dividing equation (6) 
by equation (5) 


P(n at Z/np at Zp, 


n 2 at Z 2 ) 


( n i - n 2 ) I 
(np - n) 1 (n - n 2 ) I 


_R( Zp) - R(Z)] ni ~ n [~R(Z) - R(Z 2 )] n ~ n2 
[R(Zp) - R(Z 2 )] ni ‘ n2 

The expectation value of n at Z for a bundle passing through np at Zp 
and n 2 at Z 2 is 


n l 

E(n at Z/np at Zp, n 2 at Z 2 ) = ~] nP(n at Z/np at Zp, n 2 at Z 2 ) 

n=n2 


= n 2 + (np 


R(Z) - R(Z 2 ) 

n 2 j R( Zp) - R(Ti7 


( 8 ) 


where N > n^_ > n > n 2 >0 and 0 < Zp < Z < Z 2 . Equation (8) shows that the 
average number of fibers surviving in a bundle passing through np at Zp and 
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n 2 at Z 2 is linear in R(Z). Therefore, the 
average path appears as a straight line connecting 
the points np at Zp and ng at Z 2 on a 
plot of n against R( Z ) . 


CLASSIFICATION OF PATHS 


load l, R(Z) 

Figure 2. - Average paths of bundles. 


In the calculations in this section, 
paths will be grouped in accordance with 
number of fibers surviving at Z, and the 
loads for each class of paths shall be 
function of Z or R( Z ) . It shall be 


the 

n, the 
average 
found as a 
shown that 


the average load for gaths 1 ( fig . 2 ) having more 
fibers surviving at l ^ than expected reaches a 
maximum value at 1 = 1, whereag the average load 
for paths 2 having^fewer fibers surviving than expected at Z reaches a maxi- 
mum value at Z < Z. This subdivision of paths in accordance with_the value of 
n carries through at all values of Z so that the average load L at any Z 
for bundles having a given n may be determined. The subdivision can also be 
regarded as a first, although primitive, step in the approach to the ideal 
wherein the probability and maximum supported load are determined for each 
possible path, so that an_averaging of the maximum loads over the infinite pop- 
ulation of paths yields S, the average breaking strength of the bundle. 


In accordance with equation (8), the expectation value of 
a bundle passing through n at Z is 


n at Z for 


E(n at Z/n at Z) = < 


H - (H - n) ||1 


n 


R(Z) 

R(Z) 


n > n, Z < Z 
n < n, l > l 


(9a) 

(9b) 


Figure 2 represents the situation that prevails. Every path must go through 
the initial point n = N, R( Z ) =1 and the terminal point n = 0, R(Z) = 0. 

If no other point along the path is specified, the path averaged over the 
bundles is given by curve 0 of figure 2, that is, by the straight line connecting 
the initial and terminal points. If a point is specified along the path between 
the initial and terminal points, two-segmented lines result that represent the 
average paths of bundles going through the additional point specified (curves 1 
and 2 of fig. 2). The additional point specified for^curve 1 is nj_ at Z, 
whereas the additional point for curve 2 is n2 at Z . 

When the load per fiber is Z^ the average total load supported by a 
bundle which goes through n at Z is 

L( z/n at Z) = ZE(n at z/n at Z) (10) 
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II III I 


i am n i n 


MAXIMUM LOADS ON PATHS 

^The maximum load that occurs on an average path when a particular n 
at Z is specified can be determined as follows . Let G( Z ) = ZR( Z ) ; then 
equations (l) and (3) can be written 

L(Z) = NG(Z) (la] 

G'(Z) = 0 (3a] 

The assumption is also made that the form of R^Z ) is such that L( Z ) exhibits 
one and only one peak, that peak occurring at Z. Mathematically this is 
equivalent to requiring that 

G’(Z) > 0 Z < Z (lla; 


G* ( Z ) < 0 


Z > Z 


From equations (9) and (10), 


dL( Z/n at Z ) _ 
dZ 


n - MR(Z) + (N - ri)G'(Z) Z < Z 
B(Z) 


G’(Z) 


z > z 


Equation (12a) shows that^ for values of n > NR(z), L(Z/n at l) is 
always increasing up to 1 = 1. At 1 = 1, the slope changes sign (eq. (12b)) 
and L(z/n at Z) decreases as l increases when l >1. Hence, ^when 
n > RR(Z), the maximum value of L(z/n at l) occurs at l = l or 


L( l/n at l ) r 


m(l) < n < N 


Again, by equations ^12), for values of n < N \R(l), L(Z/n at l) reaches 
a maximum when l = A£ < l, where A£ satisfies the equation 


r'(Ag) = 


HR(Z) - 
N - n 


With 5£ denoting the right side of equation (14), the following quantities 
can be expanded in powers of 5g, where the coefficients^of the powers of 5g 
are functions of the derivatives of G( Z ) evaluated at Z : 


= l + 


4 * °( 8 I) 


G(Ag) = G(l) *£-4- I % 4 + o(4) 
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(17) 


L(Ag/n at l) = nZ - 


N 

2G 2 


J n 


n(£ 3_ _ i\ 8 g 

2 \3G| g 2/ n 


0(8g) 


In equations (15) to (17) the symbol G n represents the n^* 1 derivative of 
G(Z) with respect to 7 evaluated at 1 = 1. Note that G2 < 0. 

Equations ( 13 ) and ( 17 ) show that to first order in 6£, when the paths 
are classified in accordance with n at 7 , Daniels' result holds , that is, 
the breaking strength of the bundle is only dependent on the number of survivors 
at 7 . 


Equation ( 17 ) may be employed to find a larger lower bound for S, the 
average breaking strength, than that given by equation (2): 

N 

S > P(n at 7 )L( Ag/n at l) > MR (7) 
n=0 

[mil)] 

+ 1 ^ 2 / P(fi at ^ ["^ ^ + (i ^ + * * •] (18) 

n=0 ’ 


If the variable is changed from n to y, where y is of standard measure 
and is defined by 


mji) - n 

Vnr(Z)B(Z) 


(19) 


the frequency function f(y) can be expanded in terms of Hermite polynomials 
H n (y) by the method of Kendall and Stuart (ref. 3): 


f(y) 



e-y 2/ 2 [j. + | + . . .] 


^ — - — e -y^/2 1 + — Ml) ~ 1 

L 6 VwR(7)B(z) 


(y 3 - 3y) + 


• • 


Also can be -written as a series in powers of y: 


00 

_ JsBL y V r pklkli 1 k/2 yi 
^ - M m(i) v /_J 11 ^ 

k=0 


( 20 ) 


( 21 ) 
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Substituting in equation (18) yields 


S > HtR(t) + 


GO 

1 f 1 e -y 2 /2 |i + 1 Mil 

2 J L 6 


- — (y 3 - 3 y) + . . . 

\^r(Y)b(z) 


] 


f-JL £(Xl y 2 + + M)] 

\ °2 B(TT Ug| GbJLbHTJ 


3/2 ^ 3 ^ 


+ • 


. !jdy 


>mR(z)-j-|^ri + | ;h(z)-i - + b s ; 

4G 2 B(Z) [ 3 V2jtNR(Z)B(Z) 


■] 


B*t) 


■^2nN 


■ . ^ v -w 

sUl 

U(Z). 


"|°/2 

[l + . . .]+ 0(N“b 


( 22 ) 


The variance of the breaking load based on the previous approximations is given 
as 

° 2 - ™ J)(?)2 + [i + • • •] + o(i) 


(23) 


CONCLUDING REMARKS 

A method of determining a lower bound for the average breaking strength of 
a bundle of classical fibers has been presented which involves the subdivision 
of the bundles in accordance with the number of unbroken fibers at a given 
value of load per fiber. The method also yields the approximate statistical 
distribution of the breaking strengths so that the variance of strengths may be 
calculated. The result agrees with the literature. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, September 13, 1965. 
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